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EFFECT OF COSMIC UV BACKGROUND ON STAR FORMATION IN HIGH REDSHIFT GALAXIES 
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ABSTRACT 

The effect of the cosmic UV background on star formation in high redshift galaxies is explored by means of 
high resolutions cosmological simulations. The simulations include star formation, 3D radiative transfer, and 
a highly detailed ISM model, and reach spatial resolution sufficient to resolve formation sites for molecular 
clouds. In the simulations the local radiation field in the Lyman-Werner band around star-forming molecular 
clouds dominates over the cosmic UV background by a factor of 100, similarly to the interstellar radiation 
field in the Milky Way and in a few high redshift galaxies for which measurements exist. The cosmic UV 
background, therefore, is essentially irrelevant for star formation in normal galaxies. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - stars:formation - methods: 
numerical 



1. INTRODUCTION 

Since the end of cosmic reionization, the intergalactic space 
has been filled with the accumulated ultra-violet radiation 
from the previous generations of massive stars and quasars, 
the so-called cosmic UV background. The UV background 
completely controls the ionization state of the intergalactic 
medium (IGM) that manifests itself in the numerous Ly-ar ab- 
sorpti on lines in t he spectra of distant quasars, the Ly-or forest 
(see Meiksin 2009 , for a recent review). 

In the interstellar medium of local galaxies (ISM) the sit- 
uation appears to be rever sed: the Haardt-Madau model 
for the UV background ( |Haardt & Madau| |2001| l pre- 
dicts the z = radiation field at 1000A of about 2 x 
10 3 photons / cm 2 / s/ ster/ eV, some 500 smaller than the in- 
terstellar radiation field in the Milky Way at the solar circle 
( |Draine|1978l|Mathis et al.|1983) . 

At higher redshift the situation is, however, less clear. For 
example, the UV background at z ~ 3 is expected to be 50 
to 100 times higher than at z = ( |Haardt & Madau||2001 
Faucher-Giguere et al.|2009"] l. On the other hand, estimates of 



the interstellar radiation field in high redshift gamma-ray burst 
hosts also give some 100 times higher values than the Milky 
Way field (Chen et al. 2009), but the observational constraints 
remain sparse. 

In the cosmological simulation community, the effect of the 
cos mic UV background is often included in the cooling rates 
(e.g. | C en & Ostriker 1992 ljKatz et al.|1996l|N avarro & Stein- 
metzl 1 997| |Kravtsov|2003l JCeverino & Klypin|2009j jSchaye 
et al.|2009) |Hambrick et al.||2009[ etc). However, the effect 
of the local interstellar radiation field (that dominates over the 
cosmic UV background by a larger factor at least at z — 0, and 
may be at high redshift too) has not yet been included in those 
simulations. In this paper the effect of the cosmic UV back- 
ground is critically reassessed with numerical simulations that 
both include the full 3D, time-dependent and spatially vari- 
able treatment of radiative transfer and have high enough spa- 
tial resolution to resolve the sites of molecular clouds and as- 
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sociated star formation. 

2. SIMULATIONS AND STAR FORMATION MODEL 

The physical ingredients and computational setup for the 
simulations used in this paper have been recently described 
in great detail elsewhere (Gnedin & Kravtsov 2010a. b'l. As 
a brief reminder, the simulations have been performed with 



the Adaptive Refinement Tree (ART) code (Kravtsov 



Kravtsov et al. 2002 Rudd et al. 2008[) that uses adaptive 



1999 



mesh refinement in both the gas dynamics and gravity cal- 
culations to achieve high dynamic range in spatial scale. 

The simulations include star formation and supernova en- 
richment and thermal energy feedback, as well as a highly de- 
tailed ISM model. The 3D radiative transfer of UV radiation 
from individual stellar particles i s followed self-c onsistently 
with the OTVET approximation ( [Gnedin & Abel|2001) . The 
simulations incorporate non-equilibrium chemical network 
of hydrogen and helium and non-equilibrium, metallicity- 
dependent cooling and heating rates, and a phenomenological 
model of molecular hydrogen formation on and shielding by 



cosmic dust, as well as H2 self-shielding ( Gnedin & Kravtsov 
|2010a|b) . 

Particular simulations used in this paper model a small re- 
gion including a Milky-Way progenitor galaxy and a number 
of smaller galaxies with the mass resolution of 1.3 x 10 6 M Q 
in dark matter, 2.2 x 10 5 M in baryons, and with the spatial 
resolution of 65 pc x [4/(1 + z)] (in physical units) within the 
fully refined region. 

Star formation in the simulation is oc curring in the molec- 
ular ga s o nly, using the prescription s of Krumholz & McKee 
( |20031 ) and |Krumholz & Tan| ( |2007] l. The exact formu lation of 
the star formation re cipe is shown in Equation (2) of |Gnedin] 
|& Kravtsov| ( |2010b) . 

For the purpose of this paper, two simulations are consid- 
ered that differ only by the inclusion of the cosmic UV back- 



ground from Haardt & Madau (2001 ). In the simulation with- 
out the cosmic U V background, the local radiation field pro- 
duced by nearby massive stars is still included in exactly the 
same manner as in the simulation with the UV background. 
Thus, these two simulations can be used to evaluate the par- 
ticular effect of the cosmic UV background on the properties 
of model galaxies. Because of computational expense, the 
simulations are not continued beyond z = 2. 
The simulation with the cosmic UV background is the same 
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Fig. I . — The stellar mass (top) and the star formation rate (bottom) vs the 
total mass for all fully resolved model galaxies at z = 3. Open squares show 
the run without the cosmic UV background, while filled squares show the 
fully self-consistent run, with the UV background. The star formation rates 
are averaged over 20 Myr. 
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Fig. 2. — Comparison between the proximity zones (filled diamonds) and 
radii containing 90% of the neutral (top) and molecular (bottom) gas in model 
galaxies at z = 3 (open squares). Vertical error-bars give 10% - 90% range 
for the proximity zone sizes as viewed in different directions from the center 
of each galaxy. 



one as described in Gne din & Kravtsov| ( |2010a] l. That refer- 
ence also shows good agreement of that simulation with sev- 
eral observational constraints on the properties of high red- 
shift galaxies. 

3. RESULTS 

A direct comparison of stellar masses and star formation 
rates between the galaxies in the two simulations is shown in 
Figure [TJ In order to minimize the effect of finite numerical 
resolution, only highly resolved galaxies - i.e. galaxies that 
reach the 8th level of mesh refinement (130pc spatial resolu- 
tion at z = 3) within their gaseous disks - are shown in this and 
all subsequent figures. At this spatial resolution the sub-cell 
model for H2 forma tion perform s reliably, as is demonstrated 
in Figure 13 of |Gnedin & Kravtsov| ( |20 1 0b] > . Galaxies with 
masses below M tot = 10 10 M o are not sufficiently resolved in 
these simulations - none of such galaxies achieves 8 levels of 
mesh refinement within their gaseous disks. 

While the two simulations do not produce identical results, 
the difference between the two runs is not dramatic and is fully 
consistent with the timing differences in two simulations that 
have slightly different time-steps and, hence, times of inter- 
mediate outputs. Thus, no significant effect of the cosmic UV 
background on the global properties of simulated galaxies is 
observed in these simulations. 

The origin for this conclusion becomes apparent from Fig- 
ure [2] where the sizes of the proximity zones - regions around 
model galaxies that are dominated by the local radiation from 
the galaxies themselves and not by the cosmic UV back- 
ground - are shown together with radii containing most of 
atomic and molecular gas. 

A size of a proximity zone depends, of course, on the wave- 
length of radiation, because, in general, spectral shapes of 
the local interstellar radiation field and the cosmic UV back- 
ground are different. Figure[2] therefore, shows two important 



wavelengths: the Lyman limit and A = 1000 A, in the middle 
of the Lyman-Werner band. Radiation above the Lyman limit 
ionizes neutral hydrogen, and so is important for determining 
which part of the galactic gas can cool rapidly. Radiation in 
the Lyman-Werner band destroys molecular hydrogen; only 
when gas is shielded from that radiation by cosmic dust and 
molecular hydrogen self-shielding can it become fully molec- 
ular (and, hence form stars). 

In a given frequency band, absorption of radiation emitted 
by stars in a galaxy will, generally, be different in different 
directions. Therefore, the proximity zone of a given galaxy is 
not necessarily spherical, but can have a complex shape. In 
order to quantify variatio n in the shapes of p roximity zones, 
a HEALPix tessellation (Gor ski et al.| [2005) of the celestial 
sphere (as seen from the center of a particular galaxy) is con- 
structed and the proximity zone size is measured separately 
for each pixel on the sky. Because HEALPix provides uni- 
form sampling of all possible directions, it offers a convenient 
way to quantify variations in sizes of proximity zones in dif- 
ferent wavebands. Error-bars in Figure|2]show the 10% - 90% 
range in the size of the proximity zone as a function of the an- 
gle on the sky. 

To evaluate the effect of the local radiation field on the 
gas physics, Figure [2] also displays radii that contain 90% of 
atomic (top) and molecular (bottom) gas in the model galax- 
ies. As can be seen, the proximity zone in the Lyman-Werner 
band extends over 10 times beyond the edge of the molecular 
gas. That explains why the star formation in the simulations 
is insensitive to the presence of the cosmic UV background 
- the local radiation field in the sites of star formation (i.e. 
in molecular clouds) always dominates over the cosmic back- 
ground by a factor of about 100. 

The situation appears to be reversed for atomic hydrogen. 
The escape fraction for ionizing radiation is small i n these 
simulations, consistent with observational estimates ( |Gnedin| 
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Fig. 3. — Comparison between the proximity zones and radii containing 90% 
of the neutral and molecular gas for the fully self-consistent run, with the UV 
background, at three different redshifts. Symbol and panels are similar to 
Figure|2] Three shades of gray show three different redshifts, as marked in 
the figure. 



et al. 2008 ). It is particularly small along the disk, so that little 
local ionizing radiation shines on the outer parts of H I disks; 
instead, the edges of H I disks are determined by the cosmic 
background. 

Finally, to verify the robustness of these results, Figure 
presents proximity zones and 90%-mass radii for the full sel 
consistent simulation (with the cosmic UV background) at 3 
different redshifts. While there is substantial scatter between 
individual points, the trends observed at z = 3 are reproduced 
at other redshifts as well. 

4. CONCLUSIONS 

The results of this paper can be summarized in just a few 
lines: the local radiation field in the Lyman-Werner band 
around star-forming molecular clouds dominates over the cos- 
mic UV background by a factor of 100, similarly to the inter- 
stellar radiation field in the Milky Way ( |Draine|1978| |Mathis| 



et al. 1983 ) and i n a few high redshi ft galaxies for which mea- 
surements exist ( |Chen et al.| [2009). The cosmic UV back- 
ground, therefore, is essentially irrelevant for star formation 
in normal galaxies. Only in the lowest mass dwarfs, where 
the radiation can actually photo-evaporate the ISM, may the 
cosmic background affect star formation. Such galaxies, how- 
ever, contribute very little to the total star formation history of 
the universe. 

The situation is more complex for ionizing radiation - in 
that band the cosmic UV background may play an impor- 
tant role in determining the exact locations of the edges of 
extended HI d isks, in agreement with the earlier results of 
Schaye ( 2004) . Those outside regions of galactic disks, how- 
ever, contain little molecular gas and, therefore, are inert to 
star formation. 

The simulations presented here do not continue b eyond 
z - 2, and s o can not be used for testing the results of |Ham-| 
|brick et ak] ( |2009| l, who found little effect of the cosmic UV 
background on star formation in normal galaxies at z > 2, 
but a larger effe ct at lower redshifts. Simulations of Ham- 
|brick et aT] ( |2009 1, however, did not include radiative transfer 
and, therefore, failed to account for the local radiation and its 
dominance over the cosmic UV background within the galac- 
tic proximity zones. For example, if the Milky Way galaxy 
is representative of low redshift (z < 2) normal galaxies, then 
one may expect that even at low redshifts local radiation dom- 
inates over the cosmic background within the galactic ISM, 
rendering the background irrelevant to star formation. 
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